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ABSTRACT 

Today, many extrasolar planets have been detected. Some of them exhibit properties 
quite different from the planets in our solar system and they have eluded attempts to 
explain their formation. One such case is HD f 4902Gb. It was discovered by Sato et al. 
(|2005f h A transit-determined orbital inclination results in a total mass of 114 Mq. The 
unusually small radius can be explained by a condensible element core with an inferred 
mass of 67 for the best fitting theoretical model. 

In the core accretion model, giant planets are assumed to form around a growing 
core of condensible materials. With increasing core mass, the amount of gravitation- 
ally bound envelope mass increases. This continues up to the so-called critical core 
mass - the largest core allowing a hydrostatic envelope. For larger cores, the lack of 
static solutions forces a dynamic evolution of the protoplanet in the process accreting 
large amounts of gas or ejecting the envelope. This would prevent the formation of 
HD 149026 b. 

By studying all possible hydrostatic equilibria we could show that HD 149026 b 
can remain hydrostatic up to the inferred heavy core. This is possible if it is formed 
in-situ in a relatively low-pressure nebula. This formation process is confirmed by 
fluid-dynamic calculations using the environmental conditions as determined by the 
hydrostatic models. 

We present a quantitative in-situ formation scenario for the massive core planet 
HD 149026 b. Furthermore we predict a wide range of possible core masses for close-in 
planets like HD 149026 b. This is different from migration where typical critical core 
masses should be expected. 
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1 INTRODUCTION 

At the moment, there are two competing theories for gi- 
ant planet formation. In one theory, the solar nebula frag- 
ments directl y due to gravitational instability to form a 
giant planet (jcameronl fl978l ; iDecampli fc Cameron! fl979l ) . 
In the second, the core accretion scenario, a core or proto- 
planet forms through accretion of planetesimals and when 
its mass reaches some critical value the surrounding gaseous 
envelope is supposed to become unstable and to collapse 
onto the core, in the p rocess forming the giant planet 
|Perri fc Cameron! Il974l: iMizuno et all ll97Sl) . For an in- 
depth discussion see IWuchterl et al. (2000). In this paper, 
we will assume that the nebula is gravitationally stable and 
follow the second idea. 

Soon after the d i scove ry of the first exo-planet, 51 Peg b, 
bv lMavor fc Quelozl (| 19951 ) theorists have come up with sev- 



E-mail: broeg@space.unibe.ch 



eral migration theories predicting that the planets form at 
large orbital dist ances (around 4- 5 AU) and migrate inwards. 
Both continuous l|Lin et a l. 1996) and sudden migration (the 
jumping Jupiters of Weide nschillinel 19771) h ave b een pro- 
posed. On the other hand, IWuchterll (Il996l . Il997l ) showed 
that in-situ formation could occur if sufHcent amounts of 
solids and gas are available in the planets feeding-zone. 



The planet HD 149026 b was discovered by ISato et ail 

(2005) at a distance a of only 0.042 AU. Because it was dis- 
covered by both the radial velocity and the transit method, 
its mass and density are known. It has a total mass of 
114 Mm and an unusually large density. Calculations by the 
discovery-team give the most likely core mass as 67 Mm. 



We will show that the large core of HD 149026 b cannot 
be explained by migration. It has to form in-situ to allow 
such a large subcritical core-mass. 
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2 MODELLING THE EQUILIBRIUM 
ENVELOPE STRUCTURES 

Every protoplanet in our model consists of a solid core of 
constant density (p CO rc= 5500 kgm~ 3 j]] and an envelope of 
hydrogen and helium with a helium mass fraction of 0.24. 
The composition of the envelope is assumed to be protosolar. 
Following the recommended procedure for the SCVH equa- 
tions of state, that contain now high-Z elements, their effects 
are accounted for by a somewhat enhanced He-mass frac- 
tion. Heavy elements and their condensates are fully taken 
into account in the opacities to calculate radiative transfer 
efficiency in detail. 



2.1 Model equations and assumptions 

The hydrostatic equilibrium configuration of the envelopes 
is given by the well known equations of stellar structure: 

AM 
dr 
dP 
dr 
dT 
dr 
dL 
dr 
dU 
dr 

This, together with the boundary conditions (see sec- 
tion l2.2p . determines the values for the mass M, the pressure 
P, the temperature T, the luminosity L, and the gravita- 
tional potential U as a function of the radius r. G is the 
gravitational constant and V is calculated as: 



The outer radius is given by the Hill radius: 



4Tvr 2 p(P,T), 


(1) 




(2) 






(3) 




0, 


(4) 


GM 


(5) 





V = miii s 



i (Vrad, Va 



(6) 



i.e. the adiabatic temperature gradient V s or the tempera- 
ture gradient as caused by radiative energy transport in the 
diffusion approximation, V ra d - whichever is smaller. This 
corresponds to the use of zero entropy gradient convection 
and the application of the Schwarzschild-criterion, but is 
continuously differentiable across the transition region. V s 
is directly given by the equation of state, V ra d is calculated 
as: 

„ 3 kLP 
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(see lMihalas fc Weibel-Mihaladfl999l . chapter 6 for a deriva- 
tion of this formula and explanation of all parameters), k 
is the Rosseland-mean opacity and is determined by using 
tabulated values (see sect. 12.3 . 2^) ■ 

2.2 Boundary conditions 

To solve the differential equation system we specify the fol- 
lowing boundary conditions: 

• The core radius is defined using a fixed core density p c : 



_3_M 

47T p c ' 
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where a is the planet's semi-major axis, M its mass, and M* 
the mass of the host star. 

• The luminosity L is defined as the energy libration rate 
obtained for a constant planetesimal accretion rate and the 
dissipation of the planetesimal kinetic energy at the core 
surface: 



U(r c . 



GM r 



-{U-U )M, 



dr, and 



(8) 



where Uo is arbitrary. Uo is the gravitational potential at 
rmn- 



2.3 Constituent relations 

To fully specify the differential equation system, V s , p, and 
k need to be specified. We use the following equation of state 
and opacity tables. 



2.3.1 Equation of state 

V a (-P. T), and p{P,T) are interpolated from [Saumon et al.l 
( 1995). First hydrogen and helium are interpolated indepen- 
dently, then the mixed quantity is determined. V s is calcu- 
lated via spline derivatives from the mixed e ntropy includ- 
ing th e mixing entropy term, as suggested in lSaumon et al.l 
(1995|). The helium mass-fraction is Y = 0.24. 



2.3.2 Opacity 

Rosseland-mean opacities n(p, T) are interpolated from 
a combined table: Opacities in cl ude Rosseland-mean 
dust opacities from [Pollack et all l| 19851 . lgT < 2.3), 



lAlexander fc Ferguson (1994) values in the molecular range, 
and Weiss et al. I ill990h Los Alamos high temperature opac- 
ities. 



3 UTILIZING ALL EQUILIBRIA TO 

DETERMINE THE ENVIRONMENTAL 
PROPERTIES FOR HD 149026 b 

In order to solve the system of differential equations, a range 
of parameters must be provided, namely the: 

i) core mass M co 

re , 

ii) pressure at the core P core , 

iii) mass of the host star M», 

iv) semi- major- axis of the planet a, 

v) nebula temperature T nc b, and the 

vi) planetesimal accretion rate M . 



1 a core density of 10500 gives similar results 



Parameters (i) and (ii) are our independent parameters, they 
are varied in a scale-free way, i.e. equidistant in the loga- 
rithm. 
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(iii)|(iv) and (v) are determined by the host star and 
the position of the planet. The nebula temperature can be 
calculated in thermal equilibrium with the star: 

1/4 



280 ■ 



/ a \-i/2 
V1AU7 



L 



0) 



fsee lHavashilll98ll ; lHavashi et al.lll985l ). 

The only remaining free parameter is the planetesimal 
accretion rate M. In the case of HD 149026 b we choose an 
unusually large number of M = 10 -2 M^a -1 . This value 
corresponds to particle-in-a-box planetesimal accretion for 
the density at the position of HD 149026 b in a min imum 
mass solar nebula (lHavashilll98ll ; lHavashi et al.| [l985) with 
a gravitational enhancement factor F g — 21. 

For HD 149026 b the correct values are therefore: 



0' 



• M* = 1.3 Mq, L, = 2.72 L, 

• a = 0.042AU, 

• T ncb = 1754 K, 

• M = 10" 2 M e a- 1 . 

Using these values we have calculated all hydrostatic en- 
velope structures by solving the equation system described 
in chapter [5] for a wide range of core masses and pressures. 
Using this solution manifold it is easy to determine the cor- 
rect environment that allows such a large critical core mass 
of roughly 70 M® (see figure [TJ . 

The results show that a hydrostatic equilibrium of 
gaseous envelope and solid core is indeed possible at the 
specified position, if the nebula pressure is P = 10 3 ' 6 Pa = 
3980 Pa. In that case the envelope mass is lgM onv kg _1 = 
26.45 or roughly 47 M® leading to a total of 117 M® for 
the protoplanetlj This is very close to the observed mass of 
114 M®. 

It should be noted, that this is the case only very close 
to the host star. With increasing distance, the accretion rate 
M decreases, and even more importantly the hill radius in- 
creases for constant planet mass. Both effects lead to a re- 
duced critical core mass. Therefore, such a large core is only 
allowed for very close-in planets like HD 149026 b. 



4 FLUID-DYNAMIC FORMATION OF 
HD 149026 b 

Knowing the environmental conditions, especially the values 
of the nebula pressure P ne b and accretion rate M we tried 
to reproduce the in-situ formation of HD 149026 b with a 
thorough fluid-dynamic calculation. 

The algorithm for the dyna mical calculations is de- 
scribed in IWuchterll (|l990l . Il991al lb]l For this paper we use 



a mod ified convection theory as in lWuchterl fc Tscharnuterl 
(2003). These c alculations use a different equation of state 
(not SCVH, see lWuchterill989l ) and slightly different dust 
opacities, namely for inte rstellar dust instead of the Polla ck 
et al. protosolar dust, see I Wuchterl fc Tscharnuterl |2003). 

Because the dynamical calculation uses a particle-in-a- 
box accretion scheme instead of a constant accretion rate, 



Repeating these calculations with a core density of p CO rc = 
10500 kg m — 3 and a critical core mass of 67Mgj leads to the same 
nebula pressure (P = 10 3 ' 6 Pa), a slightly larger envelope mass 



HD1 49026b: Mass and Luminosity CRFD+PiB(10) 



"E 1 250- 




(lg M env kg 



26.5, or 53Mfg) and a total mass of 120 Mq 



Age/[a]wuchterl and Tscharnuter 2003 

Figure 2. Possible evolution of HD 149026 b. The figure shows 
the time-evolution of three quantities: the core mass (dashed), the 
total mass (solid), and the luminosit y (solid). The age is defined 
as in lWuchterl fc Tscharnute j (|2003h . CRFD+PiB(10) stands for 
fully hydrodynamic simulation with Safronov number 10. The 
smooth luminosity curve is a strong hint for quasistatic evolution 
which is confirmed by looking at the data. The dynamic terms 
are negligible. 



we use a gravitational enhancement factor F g = 21 (or a 
Safronov number 9 — 10) which leads to a peak accretion 
rate of M = 10~ 2 A/® a" 1 as required. 

In spite of the slight differences regarding the con- 
stituent relations, the dynamic calculations confirm the hy- 
drostatic model. The evolution of HD 149026 b is plotted 
in figure The dynamic calculation confirms quasistatic 
evolution for the entire formation process of HD 149026 b, 
no instabilities occur. This formation scenario explains the 
high core mass and shows no dynamic accretion phase - the 
planet grows hydrostatically all the way to its final mass. 



5 SOLVING THE FEEDING-ZONE PROBLEM 

So far we have shown how HD 149026 b could have formed 
provided that there is enough material available to form the 
planet at its current position. The lack of building material 
is usually considered the strongest argument against in-situ 
formation of close-in planets. 

It is true that in a classical feeding zone, i.e. 3-4 hill 
radii (equation (0) on both sides of the orbit of the planet, 
there is not enough material available in a gravitationally 
stable disk. This problem can be circumvented by assuming 
a continuous flow of material onto the star as is the c ase for 
accretion disks. According to Ha rtmann et all ||1998() typi- 
cal T Tauri disks with an age of 1 Ma have relatively low 
accretion rates, they say: 

" The median accretion rate for T Tauri stars of age ~ 1 Ma is 
10 — 8 Mq a -1 ; the intrinsic scatter at a given age may be as large 
as 1 order of magnitude." 
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Figure 1. The figure shows the envelope mass (logM en v) as a function of core mass and pressure at the core surface (in the figure 
logMc and logPc, respectively). The logarithms are taken from the corresponding values in SI units (kg / Pa). The results are given as 
a surface in three dimensions and the surface colour is mapped from the outside pressure (logPx). The blow-up on the right-hand-side 
shows the region of interest around a core mass of 70 Mq . 

To determine the nebula conditions for a critical core mass of 70 M^, we have drawn lines of constant nebula pressure (isobars, yellow, 
in the range lgPPa - 1 = 3. .4, step 0.1) and lines of constant envelope mass (gray, in the range lgM cnv kg — 1 = 26. .27, step 0.1). The 
critical core mass for a given nebula pressure is given as the largest possible core mass for a given isobar. So the isobar that is tangentially 
touched by the line of 70Mgj determines the nebula pressure for which the critical core mass is 70Mgj. It is lgPPa - 1 = 3.6. 
Using the gray lines we can immediately give the envelope mass corresponding to this critical core mass. It is 47 giving a total mass 
of roughly 117 Mm. This is almost exactly the total mass of HD 149026 b. 



For now, we will assume a very fast formation of the 
planet of 10 a and calculate the amount of material that 
passes the orbit of the hypothetical in-situ planet. Using 
Hartmann's estimate, the material passing a close-in planet 
during its formation is therfore M = 10 5 x 10~ 8 Mq = 
333 Mq. For the envelope this is quite enough mass, but 
what about the heavy element core? For solar compo- 
sition gas we have a mass-fraction of condensible mate- 
rial of ~ 1/56 outs ide of the ice- line and « 1/240 inside 
jHavashi et al.lll985T l. Typical disk-lifetimes are ~ 1 Ma. We 
can now calculate the amount of condesible material that 
passes the protoplanet's orbit during the lifetime of the disk: 

M cond = 10 6 x 1O" 8 M /24O « 14M . (10) 

Keeping in mind that HD 149026 b is an extreme case, it is 
quite appropriate to use the upper limit given by Hartman: 



a disk accretion rate of 10 7 Mq a 1 . This provides 140 Mq 
of condensible material at the orbit of the planet. 

In the case of HD 149026 b we can go even fur- 
ther. HD 149026 's metalicity is given as [Fe/H]=0.36 (from 
ISato et al]l2005l ). Assuming that the other heavy elements 
are similarly enriched, this is an overabundance in heavy el- 
ements of lO 36 = 2.3. In total we end up with ~ 300 Mq in 
heavy elements passing the planet's orbit - This is enough 
solid material for HD 149026 b. 

There is one problem remaining, concerning the conden- 
sibles: the high temperature of the nebula. At a temperature 
of 1754 K and a pressure of 4000 P a will there be any con- 
densed specimen left? According to lDuschl et al. | (|l996l ) the 
mass fraction of silicon for this P-T-regime is 0.99-0.999. 
So the silicates are still available. For other specimen, es- 
pecially carbon, this is usually not the case. To answer the 
question of condensible mass-fraction precisely, the exact en- 
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vironmental conditions, such as pressure, temperature, and 
chemical composition at the time of HD 149026 b's forma- 
tion need to be taken into account. Without this informa- 
tion, we can only speculate. As it is very likely that cm- 
sized and larger grains spiral inward to the star, a fraction of 
the particles will never reach equilibrium before meeting the 
protoplanct. Once inside the planet, the high pressure en- 
vironment prevents evaporation. Therefore we assume that 
a fraction of the non-silicates can also be accreted onto the 
core of the protoplanet providing just enough material to 
form HD 149026 b. 



6 NOTES ON MIGRATION 

The reader might have pondered the lack of migration in 
the above discussion - this is partly intended. We wanted to 
show that while different types of migration can take place 
for different embryo masses and nebula proberties, they are 
not strictly necessary for the formation of close-in planets. 

In the case of HD 149026 b we can go even further: Our 
calculations show that to obtain such a large core, at least 
the last phase of core accretion must have occured in close 
proximity to its present location, i.e. without migration. The 
'last phase' in this case refers to the time when the core mass 
grows beyond ~ 30 Mg. Before that time, migration could 
have occured - this is irrelevant for the formation scenario 
presented here. Once the planet nears its final core mass, the 
nebula pressure is already very low. Therefore the planet's 
orbit should be stable against type-I migration. 

This scenario has some interesting consequences. The 
conventional accretion model for planet formation at larger 
separations predicts a typical core size. Two giant planets 
forming in the same nebula should be of roughly similar 
mass, i.e. giant planets orbiting the same star can be ex- 
pected to have similar core masses. This is not true for the 
in-situ formation scenario presented in this paper. In the 
case of HD 149026 b we could show that there is no dynam- 
ically triggered gas-accretion phase that sets in beyond a 
critical core mass. If this is the case for all close-in plan- 
ets, core masses will only be determined by the amount of 
available material. Therefore we expect a wide distribution 
of core masses for close-in planets in the case of in-situ for- 
mation. This property could be used to distinguish between 
in-situ formation and large-scale migration. 



7 CONCLUSIONS 

We have demonstrated a plausible formation scenario for 
HD 149026 b. In-situ formation permits the growth of the 
planet's very large core up to its present size. 

Provided that enough material is available, even very 
low densities (e.g. minimum mass solar nebula values), pro- 
duce very high planetesimal accretion rates. A safronov 
number of 10 is enough to produce accretion rates of M — 
W' 2 M a' 1 building a 70 M core in less than 10000 a. By 
analysis of all possible hydrostatic envelopes around such a 
large core we could determine the correct nebula pressure 
and demonstrate a fluid-dynamic model for the formation 
of HD 149026 b. 

In chapter [5] we analysed the mass flow in a typical 



T Tauri star accretion disk at an age of 1 Ma. We could show 
that for such a case, enough material is transported across 
the orbit of a close-in planet like HD 149026 b to provide 
sufficient material even when considering grain-evaporation 
caused by the high temperatures this close to the star. 

The planet HD 149026 b could indeed have formed in- 
situ, i.e. at its present position. In that case its entire evo- 
lution would have been quasi-static - no dynamic accretion 
takes place. Forming HD 149026 b at larger distances in this 
manner - e.g. outside the ice-line - is not possible as the 
critical core mass decreases dramatically with increasing dis- 
tance. 

Our model predicts a wide range of core masses for 
close-in planets while migrating planets should all have sim- 
ilar 'typical' critical cores. This property could be used to 
distinguish between the two formation scenarios. 



ACKNOWLEDGEMENTS 

This research was supported in part by DLR project number 
50-OW-0501. 



REFERENCES 

Alexander D. R., Ferguson J. W., 1994, ApJ, 437, 879 
Cameron A. G. W., 1978, Moon and Planets, 18, 5 
Decampli W. M., Cameron A. G. W., 1979, Icarus, 38, 367 
Duschl W. J., Gail H.-P., Tscharnuter W. M., 1996, A&A, 
312, 624 

Hartmann L., Calvet N., Gullbring E., D'Alessio P., 1998, 
ApJ, 495, 385 

Hayashi C, 1981, Progress of Theoretical Physics Supple- 
ment, 70, 35 

Hayashi C, Nakazawa K., Nakagawa Y., 1985, in Proto- 
stars and Planets II, pp. 1100-1153 

Lin D. N. C, Bodenheimer P., Richardson D. C, 1996, 
Nature, 380, 606 

Mayor M., Queloz D., 1995, Nature, 378, 355 

Mihalas D., Weibel-Mihalas B., 1999, Foundation of Radi- 
ation Hydrodynamics. Dover Books 

Mizuno H., Nakazawa K., Hayashi C, 1978, Progress of 
Theoretical Physics, 60, 699 

Perri F., Cameron A. G. W., 1974, Icarus, 22, 416 

Pollack J. B., McKay C. P., Christofferson B. M., 1985, 
Icarus, 64, 471 

Sato B., Fischer D. A., Henry G. W., Laughlin G., Butler 
R. P., Marcy G. W., Vogt S. S., Bodenheimer P., Ida S., 
Toyota E., Wolf A., Valenti J. A., Boyd L. J., Johnson 
J. A., Wright J. T., Ammons M., Robinson S., Strader J., 
McCarthy C, Tah K. L., Minniti D., 2005, ApJ, 633, 465 

Saumon D., Chabrier G., van Horn H. M., 1995, ApJS, 99, 
713 

Weidenschilling S. J., 1977, Ap&SS, 51, 153 

Weiss A., Keady J. J., Magee N. H., 1990, Atomic Data 

and Nuclear Data Tables, 45, 209 
Wuchterl G., 1989, PhD thesis, Univ. Wien 
— , 1990, A&A, 238, 83 
— , 1991a, Icarus, 91, 39 
— , 1991b, Icarus, 91, 53 



© 2002 RAS, MNRAS 000, [Lj{6] 



6 C. Broeg and G. Wuchterl 



— , 1996, Bulletin ol the American Astronomical Society, 
28, 1108 

— , 1997, in Science with the VLT Interferometer, Paresce 

F., ed., Springer, Berlin, pp. 64 — h 
Wuchterl G., Guillot T., Lissauer J. J., 2000, Protostars 

and Planets IV, 1081 
Wuchterl G., Tscharnuter W. M., 2003, A&A, 398, 1081 



© 2002 RAS, MNRAS 000,(21] 



